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Plasma  excitations  have  been  of  considerable  interest  in  recent  studies  of  semiconduc- 
tor quantum  wells  and  quasi-two-dimensional  (2D)  conduction  channels  in  FETs  and 
HEMTs.  The  intrasubband  plasmons  that  are  associated  with  a single  electron  subband 
can  be  studied  as  collective  excitations  in  a two-dimensional  electron  fluid.  When  the 
electron-electron  collision  time  is  much  smaller  then  the  collision  times  with  impurities 
and  phonons  the  hydrodynamic  model  should  be  applicable  to  the  carriers  in  MOS- 
FET  and  HEMT  channels.  Various  non-linear  wave  phenomena  were  predicted  for  the 
electron  fluid  [1],  Recently,  Shur  and  Dyakonov  [2]  analyzed  new  effects  related  to 
plasma  oscillations  and  proposed  novel  electronic  devices  operating  in  terahertz  fre- 
quency range.  The  velocity  of  linear  plasma  waves  vq  is  determined  by  the  gate  voltage 
swing  Uq,  sq  = ( eUo/m )’/2  where  m is  the  electron  effective  mass,  e is  electron  charge. 
In  the  case  of  GaAs  HEMT  for  the  gate  lengths  from  1 /im  to  0.1  /im  at  U$  = 1 V 
the  fundamental  plasma  resonance  frequency  u>o  varies  from  0.5  to  5 THz  [2].  It  was 
shown  in  ref  [1,  2]  that  under  asymmetric  boundary  conditions  on  the  source  and  drain 
contacts  the  linear  plasmons  become  unstable  and  new  non-linear  effects  in  the  device 
response  appear. 

In  this  work  we  consider  a non-linear  device  response  to  harmonic  signals  applied 
as  AC  electric  potentials  at  the  source  and  drain  contacts.  We  consider  one  dimensional 
fluctuation  of  2D  density  along  the  channel  (x  axis).  Let  «o  be  an  equilibrium  2D 
density  of  electrons  in  the  channel,  determined  by  the  gate  voltages.  We  write  the  time 
and  space  dependent  2D  density  as  n(x,  t)  = n0  + Sn(x,  t).  The  basic  equations  are  [2] 


dn  d(nv ) 

+ =0 
dt  dx 

(1) 

dv  dv  e dU  v 

dt  dx  m dx  r 

(2) 

where  dU /dx  is  the  electric  field  component  along  the  channel,  v(x,  t)  is  local  velocity 
of  the  electron  fluid,  r is  an  electron  collision  time  with  phonons  and  impurities.  If 
the  distance  d from  the  channel  to  the  gate  is  small  compared  to  the  channel  length  L 
and  the  edge  effects  are  neglected  the  Poisson  equation  which  relates  U and  n can  be 
treated  in  gradual  channel  approximation  resulting  in  a linear  local  relation: 

C 

n(x,t)  = —U(x,t)  (3) 
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Fig  1.  Dimensionless  source  current  density. 


where  C is  the  gate  capacitance  per  unit  area.  The  deviations  from  this  approximation 
in  higher  orders  of  d/L  lead  to  the  non-linear  dispersion  effects  in  plasma  waves  [3] 
and  are  not  included  here. 

In  our  numerical  solution  we  introduce  dimensionless  variables:  length  x/L,  time 
t/T  where  T = L/sq,  density  n/n0,  velocity  v/sq,  frequency  u>/u>q.  We  also  introduce 
a dimensionless  friction  coefficient  7 = L/sqt.  The  r itself  can  be  determined  from 
a known  mobility  in  the  channel.  We  take  boundary  conditions  corresponding  to  the 
application  of  harmonic  signals  at  the  source  and  the  drain.  We  want  to  see  how  the 
high  frequency  plasma  oscillations  lead  to  a non-linear  response  when  applied  signals 
are  of  much  lower  frequency.  Prom  equation  (3)  the  boundary  conditions  for  applied 
electric  potentials  can  be  stated  in  terms  of  density  variations:  5n(  0,  t)/»o  = As  cos  (u>st), 
5n(  1,  f)/«o  = Ad  cos(tOdt).  From  the  numerical  solution  of  equations  (l)-(3)  we  find 
the  source  (or  drain)  current  density  as  /,(/)  = n( 0,  f)v(0,  t).  It  is  shown  in  Fig.  1 as  a 
function  of  time  for  ^4,  =Aj  = 0.1,  u>s/u>o  = 0.01,  cu^/tuo  = 0.03,  g = 0.001.  We  see 
that  after  a short  transition  a periodic  regime  is  attained.  In  the  Fourier  analysis  of  j(t) 
many  harmonics  besides  those  of  the  applied  signals  are  generated.  The  form  of  /(/) 
is  found  to  be  determined  by  the  non-linear  plasma  waves  in  the  channel.  In  particular 
during  the  time  intervals  corresponding  to  large  drops  and  raises  of  the  source  current 
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a shock  wave  propagates  in  the  channel,  reflecting  a few  times  between  the  contacts. 
The  corresponding  density  profile  is  shown  in  Fig.  2. 

We  find  that  the  non-linear  wave  propagation  leads  to  the  interesting  non-linear 
response  in  the  terminal  current  even  though  the  applied  signal  is  of  much  lower  fre- 
quency then  the  fundamental  plasma  frequency.  Similar  response  is  found  when  the 
frequencies  of  applied  signals  are  equal  but  their  amplitudes  are  different  or  if  there  is  a 
phase  shift  between  applied  signals.  This  may  allow  for  the  experimental  demonstration 
of  the  non-linear  plasma  response  and  some  novel  device  applications. 

Acknowledgment:  G.S.  was  supported  by  the  U.S.  Army  European  Research  Office 
(London),  contract  No.  N68171-97-M-5548. 

References 

[1]  M.I.  Dyakonov  and  M.S.  Shur,  Phys.  Rev.  Lett.  71  2465  (1993). 

[2]  M.I.  Dyakonov  and  M.S.  Shur,  IEEE  Transactions  on  Electron  Devices  43  380  (1996). 

[3]  S.  Rudin  and  G.  Samsonidze,  in  Proceedings  of  ISDRS-97,  Charlottesville,  Virginia,  p.  583 
(1997). 


